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ABSTRACT 

We have obtained deep optical spectroscopic data of the highest-redshift cluster can- 



didate (z ~ 1.4, CVB13) selected by Van Breukelen et al. (2006) in a photometric 
optical/infrared catalogue of the Subaru XMM-Newton Deep Field. The data, which 
comprise 104 targeted galaxies, were taken with the DEep Imaging Multi-Object Spec- 
trograph (DEIMOS) on the Keck 2 telescope and yielded 31 secure redshifts in the 
range 1.25 < z < 1.54 within a 7' x 4' field centred on CVB13. Instead of one mas- 
sive cluster at z = 1.4, we find evidence for three projected structures at z = 1.40, 
z = 1.45, and z — 1.48. The most statistically robust of these structures, at z = 1.454, 
. has six spectroscopically confirmed galaxies. Its total mass is estimated at 10 14 M Q 

and it may therefore be termed a cluster. There is an X-ray source at the cluster po- 
sition which is marginally spatially resolved but whose X-ray spectrum is too hard to 
be thermal cluster emission. Its origin could be the summed X-ray emission from ac- 
tive galaxies in, and projected onto, the cluster. Serendipitously we have discovered a 
cluster at z — 1.28 with a mass of 10 14 M Q at the same position on the sky, compris- 
ing six spectroscopically confirmed cluster galaxies and at least one additional radio 
source. The selection of CVB13 for the cluster catalogue was evidently aided by the 
. superposition of other, presumably lower-mass, structures, whereas the single cluster 

at z = 1.28 contained too few galaxies to be isolated by the same algorithm. Given the 
complicated nature of such structures, caution must be employed when measuring the 
mass function of putative high-redshift clusters with photometric techniques alone. 

Key words: galaxies: clusters: general - galaxies: high-redshift - radio continuum: 
general - X-rays: galaxies - X-rays: galaxies: clusters 
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1 INTRODUCTION 

Clusters of galaxies are important probes for both cos- 
mology and astrophysics. They are the most massive viri- 
alised structures in the Universe, and they are therefore 
excellent targets to study large-scale structure formation. 
The cluster mass function directly links to the normalisa- 
tion of the power spectrum of density perturbations, o~$, 
which is the RMS in the density on a scale of 8/1 -1 Mpc 
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(h = Ho/WO kms _1 Mpc _1 ). The evolution of the mass 
function with redshift reflects the growth function of the 
Universe, which is determined by Qm (the matter density 
parameter), Qa (the dark energy density parameter), and 
w (the dark energy equation-of-state parameter). Further, 
clusters are of great interest to astrophysics because they can 
be approximated by a 'closed-box' environment. This makes 
them ideal laboratories to study the interaction of galaxies 
and the intergalactic medium, which is thought to play a 
large role in galaxy formation and evolution. An important 
remaining question is why so little baryonic matter in clus- 
ters forms stars (e.g. Cole 1991; White & Frenk 1991). Feed- 
back processes coupling baryonic matter from galaxies to the 
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intergalactic medium are likely to counter further conden- 
sation of the intergalactic gas and halt star-formation. Two 
examples of these processes are galactic winds due to mul- 
tiple supernovae in star-burst galaxies (e.g. Heckman 2002) 
and outflow from radio-loud AGN (e.g. Scannapieco & Oh 
2004; Fabian, Celotti & Erlund 2006a). 

Cluster searches have been carried out for decades. The 
first extensive cluster catalogue was created by Abell (Abell 
1958; Abell et al. 1989). Unfortunately Abell only had single- 
filter photographic plates at his disposal and therefore his 
work was severely complicated by projection effects. Glad- 
ders et al. (1998) improved upon this situation by showing 
that a redshift estimate can be determined by using only two 
filters and targeting the bright, red elliptical galaxy popu- 
lation in clusters which form the red sequence (e.g. Bower, 
Lucey & Ellis 1992). A particularly important advance in 
optical cluster detection has come from the advent of large 
arrays of CCD detectors allowing large surveys to be carried 
out efficiently, such as the relatively shallow (z < 0.4) but 
very wide-field Sloan Digital Sky Survey (SDSS) (e.g Goto 
et al. 2002; Kim et al. 2002; Miller et al. 2005). There have 
been numerous smaller-area surveys to higher redshifts, for 
example the Palomar Distant Cluster Survey (Postman et al. 
1996), the ESO Imaging Survey (Lobo et al. 2000), and the 
Red Sequence Cluster Survey (Gladders et al. 2005). Mov- 
ing to even higher redshifts, proto-clusters have been found 
by focussing on fields around quasars or radio galaxies (e.g. 
Venemans et al. 2002) as a large fraction of these objects 
have been shown to reside in clusters (see Pascarelle et al. 
1996; Rottgering et al. 1996). However, selecting clusters 
in blank sky surveys with optical photometry at redshifts 
z <; 1 was impractical for a long time because of the shifting 
of the 4000 A break - characteristic for the red, passively 
evolving ellipticals predominantly found in clusters - out of 
the /-band to longer wavelength bands. 

The launch of the ROSAT satellite in 1990 greatly ad- 
vanced the study of clusters in the X-ray regime, enabling 
the discovery of hundreds of clusters up to z ~ 1. Later 
Chandra and XMM-Newton provided the possibility to ob- 
serve even deeper owing to their unprecedented sensitivity 
and angular resolution. X-ray cluster surveys now extend 
well above z ~ 1 (e.g. Stanford et al. 2001; Rosati et al. 2004; 
Bremer et al. 2006) . Stanford et al. (2006) currently hold the 
record of the most-distant spectroscopically-confirmed clus- 
ter at z = 1.451. This cluster was initially detected as an 
extended source in XMM X-ray data; spectroscopic follow- 
up revealed five cluster galaxies within a l'-diameter. The 
temperature of the intra-cluster medium (ICM) was shown 
to be ~ 7 keV, which is the highest detected in any cluster 
at z > 1, implying a relatively massive cluster for such a 
high redshift. 

Another way to observe the ICM is through the 
Sunyaev-Zel'dovich effect (S-Z effect, Sunyaev & Zel'dovich 
1970, 1972). This effect is visible as a distortion of the Cos- 
mic Microwave Background (CMB): as the CMB-photons 
travel through the ICM they are subjected to inverse Comp- 
ton scattering, which shifts the energies of a small fraction 
of the photons slightly upwards. In the last decade many de- 
tections have been made of the S-Z effect (e.g. Birkinshaw, 
1999; Carlstrom et al. 2000, Jones et al. 2005) and com- 
prehensive S-Z surveys are underway (Kneissl et al. 2001). 



Cotter et al. (2002) used pointed S-Z observations to detect 
hot gas associated with az~l cluster of radio sources. 

Recent development of wide- field infrared cameras, such 
as the Wide-Field CAMera (WFCAM) on the United King- 
dom InfraRed Telescope (UKIRT), have brought photomet- 
ric cluster surveys back into vogue by enabling deep, large- 
scale infrared surveys and pushing the limit of photomet- 
ric cluster selection to significantly higher redshifts. The 
UKIRT Infrared Deep Sky Survey (UKIDSS, Lawrence et 
al. 2006) is a suite of both wide and deep infrared imag- 
ing surveys using WFCAM. The Ultra Deep Survey is the 
deepest of these, and covers 0.8 deg 2 in the Subaru XMM- 
Newton Deep Field (SXDF, Sekiguchi et al. [in prep.]) with 
a planned limiting magnitude of Ky cg!1 = 23. Deep BVRi' z' 
imaging data from the Subaru Telescope are also available 
on the SXDF (Furusawa et al. [in prep.]), as is X-ray data 
by XMM-Newton (Watson et al. 2004), radio data from the 
Very Large Array (Simpson et al. 2006, Ivison et al. 2007), 
and Spitzer infrared data ranging from 3.6 to 24 micron 
(Lonsdale et al. 2005). 

We previously exploited these exceptional multi- 
wavelength datasets to execute a high-redshift cluster sur- 
vey, using photometric redshifts, as reported in Van Breuke- 
len et al. 2006 (VB06). In VB06 we isolated 13 cluster candi- 
dates at 0.6 < z < 1.4 in 0.5 deg 2 of the SXDF. The highest- 
redshift cluster candidate, CVB13 (RA = 02 h 18 m 10.5 s , 
Dec = -05°01'05" [J2000]), was estimated to lie at z = 
1.39 ±0.07. 

In this paper we present the spectroscopic follow-up of 
the putative CVB13 cluster with the DEep Imaging Multi- 
Object Spectrograph (DEIMOS, Cowley et al. 1997) on the 
10-metre Keck 2 telescope. Section 2 describes the obser- 
vations and data reduction. The analysis and discussion on 
the nature of the system are given in Section 3. Section 4 
contains concluding remarks. Throughout this paper we use 
H = 70kms~ 1 Mpc~ 1 (h = 0.7), Q. M = 0.3, and Q A = 0.7 
and a ACDM power spectrum with as = 0.75; all sky coor- 
dinates are quoted in equinox J2000. 



2 THE DATA 

2.1 Selection of the targets 

We selected four samples of target galaxies to create a Multi- 
Object Spectroscopy (MOS) mask. Each of these samples 
comprised groups of different priorities. The samples are 
described below; the masks were created by allocating the 
maximum number of slits to sample 1, priority 1, and mov- 
ing down through all the priorities within the sample before 
allocating slits to the next sample. 

Sample 1: This sample included all galaxies that were 
determined to be a possible cluster member of CVB13 by the 
algorithm described in VB06. The algorithm presented in 
this paper used two methods to detect clusters: Voronoi Tes- 
selations and Friends-of-Friends. The galaxies that were as- 
signed to the cluster by both methods composed the highest- 
priority target group within this sample. The second prior- 
ity group comprised the cluster members as selected only 
by Friends-of-Friends, followed by the galaxies found only 
by Voronoi Tesselations, as this is a less reliable method to 
determine cluster membership (see VB06). The target list 
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Figure 1. Histogram of the photometric redshifts of galaxies in 
the field of CVB13 (02 h 17 m 54.6 s < RA < 02 h 18 m 26.4 s and 
-05°03'05" < Dec < -04°59'05" in the redshift range 1.2 < z < 
1.6. At z = 1.4 the projected region of interest has proper size 
3.6 X 2.0 Mpc 2 . The photometric redshifts were calculated using 
Bruzual & Chariot (2003) modelled spectral energy distributions 
from stellar population synthesis (see VB06 for a discussion of the 
photometric redshifts). The solid line is the sum of the redshift 
probability functions (normalised to the peak of the histogram) of 
the galaxies with a photometric redshift within the shown redshift 
range. 

of this sample consisted of a total of 125 objects, of which 
18 objects of priority 1, 3 of priority 2, and 104 of priority 3. 
Slits could be placed on 61 of these 125 targets, divided into 
priorities 1: 10 galaxies, 2: galaxies, and 3: 51 galaxies. 

Sample 2: By rotating the DEIMOS instrument by 90 
degrees we were able to utilise the long axis of the field-of- 
view to target another candidate cluster from VB06 within 
the same observations. This was cluster 6 (CVB6) at z = 

0. 76 ± 0.12, RA = 02 /l 18 m 32.7 s , Dec = -05°01'04". The 
target galaxies for this sample were selected and prioritised 
exactly as for sample 1. There were 152 targets in total, 
divided into priorities 1: 30, 2: 1, 3: 121. Of these, 64 sources 
could be observed: 13 of priority 1, 1 of priority 2, and 50 of 
priority 3. 

Sample 3: The third sample consisted of various active 
galaxies is the field, detected by their X-ray activity, radio 
emission, or excess flux in the Spitzer 24-/im band. We found 
48 radio sources in the field, which made up priority group 

1. The second priority group included 31 X-ray sources. The 
last priority group consisted of 12 objects with a detection 
in the 24- fim band and a colour of J — K — 2. The target list 
comprised 91 sources, of which 26 were observed (15 radio 
sources, 8 X-ray sources, 3 24- /im sources). 

Sample 4: The final sample was designed to fill up all 
remaining slits in the two masks. It included all galaxies in 
the field with a photometric redshift of 1.25 < z v hot < 1.54. 
Fig. [T] is a histogram of the photometric redshifts of these 
objects; overplotted is the sum of their redshift probabil- 
ity distribution functions. The chosen redshift range corre- 
sponds to the ± 2-sigma range on the photometric redshift 
of the cluster candidate CVB13. This error on the cluster 
redshift reflects the combined redshift probability distribu- 
tion functions of the cluster candidate members as given in 
VB06. This sample consisted of 1139 objects of which 43 
were observed. 



2.2 Observations and data reduction 

Optical spectroscopy was undertaken on the 21st of Decem- 
ber 2006 using the DEIMOS instrument on the Keck 2 tele- 
scope. The telescope was pointed at RA = 02 ,1 18 m 20.00 s , 
Dec = -05°0l'04.9" and the field-of-view of the instrument, 
with a position angle of 90°, was 16.7' x 5' in RA x Dec. 
The two MOS masks comprised one-arcsecond wide slits. 
The 600 line mm -1 grating was used in conjunction with 
the GG495 order-blocking filter. A long-slit observation of 
Feige 110 was taken for spectrophotometric calibration. Four 
1800-s integrations were taken for each mask, and short- 
exposure images were taken through the mask after each 
spectroscopic integration to ensure mask alignment. Seeing, 
as measured on the alignment images, varied between 0.6" 
and 1.0". Conditions were photometric throughout. 

The two-dimensional spectra were reduced with the 
DEEP2 DEIMOS Data PipelineQ- 

The pipeline processes 
the flats and arcs to determine the position of the slitlets 
on the CCD and to find the two-dimensional wavelength 
solution for each slitlet. Subsequently the spectra are flat- 
field corrected and the curvature of the spectra in the spatial 
direction is rectified. The separate science exposures are sky- 
subtracted and finally combined into one mean, cosmic-ray 
rejected, science frame. 

We extracted the one-dimensional spectra with a box- 
car extraction routine, using a 0.8"-width aperture. To en- 
able flux calibration, the standard star longslit exposure was 
reduced with the IRAF package 'twodspec'. The sensitivity 
function was determined using the IRAF task 'sensfunc' and 
the one-dimensional science spectra were calibrated with the 
task 'calibrate'. We note that for the average seeing of 0.8" 
we missed ~ 13 per cent of a typical galaxy's flux within our 
aperture. 

2.3 The spectroscopic sample 

Overall the two MOS masks targeted 194 objects. Of these, 
we were able to assign redshifts to 139 galaxies (72 per cent), 
which were divided over the samples as follows: 38 redshifts 
in sample 1, 57 redshifts in sample 2, 23 redshifts in sam- 
ple 3, and 21 redshifts in sample 4. The remaining 28 per 
cent of the observed galaxies had continua too faint to al- 
low a redshift determination from absorption features and 
showed no emission lines. Fig. [2] shows the positions of the 
targets within the DEIMOS field-of-view where each sample 
is depicted in a different colour. 

In this paper we focus solely on cluster candidate 
CVB13. We therefore imposed two selection criteria to cre- 
ate our final spectroscopic sample: i) the objects had to be 
within the 7' x 4' field centred on RA = 02 h 18 m 10.5 a and 
Dec = -05°01'05" (see Fig. [Q caption for RA and Dec lim- 
its of the sample); ii) the redshift range was restricted to 
1.25 < z < 1.54, which was the original 2-sigma range of 
cluster candidate redshift (see Section 2.1). There were 30 
objects in our DEIMOS data that satisfied these criteria, 
of which 25 were from target-sample 1 and 5 from target- 
sample 4. These objects are shown in Fig.[2]with red circles. 
We added one galaxy to our final sample that was observed 

1 http:/ /astro. berkeley.edu/~cooper/deep/spec2d/ 
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Figure 2. Spectroscopy targets overlaid on the i'-band image (Furusawa et al. [in prep.]) of the DEIMOS field-of-view. Sample 1 is 
depicted in green symbols, sample 2 in blue, sample 3 in purple, and sample 4 in yellow (see Section 2.1). The small squares symbolise 
all the objects in the target list, the plusses are the objects that were actually observed, and the diamonds are all the targets for which a 
redshift could be determined. The red circles are the objects which are included in the final spectroscopic sample for this paper: they fall 
within the depicted box of 7' X 4' centred on RA = 02 h T8 m 10.5 s , Dec = -05°01'05" and have redshifts of 1.25 < z < 1.54 (see Section 
2.3). 



with the Gemini Multi-Object Spectrograph (CMOS, Hook 
et al. 2004) on the Gemini North Telescope in Hawaii and 
satisfied the same selection criteria. The respective sample 
selection and data description can be found in Van Breuke- 
len et al. (in prep.). This brought the total of the final spec- 
troscopic sample to 31 objects. 



3 ANALYSIS AND DISCUSSION 

3.1 Properties of the line emitters 

Due to the high redshift of our spectroscopic sample of 31 
galaxies, the continua of the spectra are weak, rendering red- 
shift determination by absorption features alone unfeasible. 
All redshifts in the final sample were therefore calculated 
from the [On] emission line at 3727 A. The identification of 
the [O n] line is reliable because in some cases the doublet 
was resolved, in others there were supporting absorption fea- 
tures, or there were no other emission lines visible within the 
spectral range that would be expected if the [O n] line was 
actually another emission line at a different redshift. The 
one- and two-dimensional spectra of the detected emission 
lines are shown in Fig. 3, where the identified line features 
are labelled and noisy sky regions are shaded. The three- 
dimensional distribution of the galaxies within our target 
area is depicted in Fig. U in comoving coordinates (X cor- 
responds to RA, Y to Dec and z to redshift). The galaxies 
are colour-coded to redshift, from blue at z — 1.25 to red at 
z = 1.53. 

To determine the exact redshift, the line flux and the 
equivalent width of the [O n] emitters, we fit a line profile to 
the [On] emission line in each object. As the [On] 3727 line 
is a doublet, the fitted line profile consists of two Gaussians 
at X reat = 3726.1 A and A reat = 3728.8 A. The Full Width 
Half Maximum (FWHM) of each Gaussian is assumed to be 
equal, and is a free parameter of the fitted function. The 
other parameters are the redshift, the continuum level, and 



the ratio of fluxes of the two lines. The flux ratio is con- 
strained to be within the limits for high and low density (Os- 
terbrock 1967), corresponding to 0.35 ^ J3729/F3726 ^ 1.5, 
where F is the single line flux. Random errors on the flux 
and redshift measurement are determined by fitting mod- 
elled lines with the same procedure. 

We calculate the Star-Formation Rate (SFR) of the 
galaxies using the equation given by Kennicutt (1998): 

SFR = (1.4 ± 0.4) x 10" 34 L [o irj , (1) 

where the SFR is in Mq yr _1 and iroiij is the luminosity of 
the [0 11] 3727 emission line in Watts. This equation assumes 
L[oil]/Lua = 0.57 and a Salpeter initial mass function with 
stellar masses between 0.1 and 100 Mq. Extinction or to- 
tal obscuration is not corrected for and therefore these esti- 
mates represent lower limits on the SFR. The properties of 
the complete sample can be found in Table [T] In this table 
the redshift has been converted to heliocentric redshift (see 
Danese, De Zotti & Di Tullio 1980). 



3.2 Sample completeness 

To estimate the completeness of our spectroscopic sample, 
we have to take three factors into account: 

(i) The original target list was compiled from the 
UKIDSS-UDS Early Data Release (EDR) catalogue, which 
has a limiting if-band magnitude of ifvcga = 20.6. Cira- 
suolo et al. (2006) present a if-band luminosity function of 
galaxies in the UKIDSS-UDS in different redshift bins. We 
can calculate the total expected number of galaxies using 
their Schechter function for redshift range 1.25 < z < 1.50, 
characterised by 0* = 2.6 ± 0.3 x 10~ 3 Mpc" 3 , M£ Voga = 
-24.93 ± 0.16, and a = -0.92 ± 0.18. In our field of 7' x 4' 
and 1.25 < z < 1.54 we would expect 391 galaxies in total; 
when we impose the magnitude limit of Ky asa = 20.6 we 
would only expect to detect 99 objects, or 25 per cent. 
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Figure 3. Onc-dimensional and two-dimensional DEIMOS spectra of the galaxies in the field of CVB13. Identified emission and 
absorption features are labelled. The shaded regions denote noisy parts of the spectra due to sky lines. 
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Figure 3 cont'd. One-dimensional and two-dimensional DEIMOS spectra of the galaxies in the field of CVB13. Identified emission 
and absorption features are labelled. The shaded regions denote noisy parts of the spectra due to sky lines. 
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Figure 3 cont'd. One-dimensional and two-dimensional DEIMOS spectra of the galaxies in the field of CVB13. Identified emission 
and absorption features are labelled. The shaded regions denote noisy parts of the spectra due to sky lines. 



8 C. van Breukelen et al. 



GVHI 



ttl - 1.4EHE 






■■ t J i 






L t f i 




sH L 


'-1*1 

1 1 f 1 






Lll 


Llf 



1 
i 



6000 9 LOO 



6K» 9304 



6100 9500 



C YE 13.3? al t 




9W0 9304 9100 



6000 9 LOO 

*L i - 3,5019 



9500 




-e«nr" L 



6000 9 LOO 



6HW 9300 
VavalangLh (ii) 



9100 



95O0 



-a*io"" L 



CVB L3.3S ItZ - ].4&l7 

r 




■. j> -r_>- -. ■_ ;■. T J r — > . : : .. ™^\f.T- 1 \ i£ ^VJ" -1. 



9000 9100 WOO 9300 9400 

■ta»*l*flcLh (A? 



0500 



*.Q«i(r ,L 

B.OBlflr ,L 
-2.0*]D JI 




90O0 9100 



6800 9300 



9400 9500 



CVB13.30 aL r = l.LUOl 



i.E«mr" 

B.oicur* 1 
o.9*nr* 1 

B.aKlfflr ,L 


-S.0*]ft JI 




6000 9100 



6200 63O0 
T*v*l*ntLh (A) 



6400 9»0 



Figure 3 cont'd. One-dimensional and two-dimensional DEIMOS spectra of the galaxies in the field of CVB13. Identified emission 
and absorption features are labelled. The shaded regions denote noisy parts of the spectra due to sky lines. 



(ii) The nature of a MOS mask makes complete sam- 
pling of the target list impossible. Out of our original target 
list 177 objects satisfied our selection criteria (see Section 
2.1), using photometric redshifts for the redshift criterion of 
1.25 < z < 1.54. Of these objects, 75 were observed with the 
MOS masks. This means our spectroscopic data sampled 42 
per cent of the galaxies we intended to target. 

(iii) All the objects in our sample are [O n] emitters. 
We only detect the [On] lines with a line flux > 1.0 x 
10 34 Wm" 2 . Furthermore, sky lines raise the flux limit 
of the [O n] line that we can detect in specific parts of 
the spectra. However, this affects lines with different line 
widths in different ways. We determine our detection com- 
pleteness of [O n] lines with a flux greater than the flux 
limit of 1.0 x 10 34 Wm -2 by modelling a two-dimensional 



[O n]3727 line. The modelled lines have a rest-frame FWHM 
varying between 1.3 and 7.6 A and fluxes from 10 -20 to 
10~ 19 Wm -2 , corresponding to the range of observed rest- 
frame FWHMs and line fluxes. We add the modelled lines at 
redshifts of 1.25 z ^ 1.54 to a two-dimensional sky frame 
and try to recover the lines with the same procedure as used 
for the real spectra. Assuming a uniform distribution over 
the FWHM range and a Schechter distribution in flux, we 
calculate the fraction of missed sources at each redshift bin 
to find the detection completeness of [O n] lines with flux 
> 1.0 x 10 34 Wm -2 . The result ranges from 22 per cent at 
the locations of the strongest sky lines up to 97 per cent 
at the regions of the spectra with the least sky background 
noise. 
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ID RA Dec. z spcc F [OII] Lr OII] EW SFR i' L K 

[h m s] [° ' "] [10- 20 W m- 2 ] [10 34 W] [A] [M yr" 1 ] [L* K ] 
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Table 1. Properties of the [On] emitters. Column 1 states the galaxy ID; 24B was discovered serendipitously in the slit of 24A and docs 
not meet the A'-limit of the parent sample. The RA and Dec given in columns 2 and 3 are in equinox J2000. The (heliocentric) redshift 
and line flux of the emitters (columns 4 and 5) are derived from a double Gaussian fit to the [O II] 3727 A line profile. The line luminosity 
is shown in column 6, and column 7 is the rest-frame equivalent width. Column 8 shows the star-formation rate derived from Lron] 
(following Kennicutt 1998) and column 9 is the i'-band magnitude (Vega). Column 10 shows the if -band luminosity in units of L* at 
the galaxy's redshift, calculated with = —24.34 (Lin et al. 2004) at z = and assuming passive evolution (following prescriptions of 
Bruzual & Chariot [2003]) with a formation redshift of Zf orm = 4. The error on the flux measurement is ~ 1 X 10~ 21 Wm — 2 , although 
in regions dominated by sky lines (1.25 $ 2 ^ 1.28, 1.37 Js z ^ 1.41, and z ?J 1.5) the error is on average four times higher. The error on 
the redshift measurement is estimated at 5 X 10 -5 . The error on the luminosity is the propagated flux error. The error quoted for the 
SFR is the error derived from Eq.[T]and the error on the luminosity added in quadrature. *Observed with CMOS (see Van Breukelen et 
al. [in prep.]). 



Fig. [5] {left) shows the redshift distribution of the 31 
objects in our sample. Overplotted on the spectroscopic- 
redshift histogram is the expected number of [O n] emitters 
per redshift bin in our field, calculated with the observed 
Schechter luminosity function for [On] emitters at z = 1.47 
from Ly et al. (2006): log 10 (<^7Mpc" 3 ) = -1.97 ± 0.06; 
log 10 (L*/W) = 35.60 ± 0.05; a = -0.78 ± 0.13. The ex- 
pected number of [O n] emitters is corrected by multiplying 
by the wavelength-dependent product of the three complete- 
ness factors, namely 0.25 x 0.42 x (0.22 - 0.97). Note that 
this implicitly assumes no correlation between [O n] flux and 
observed J\~-band magnitude. However, if for example the 
rest-frame equivalent width of the [O n] line is constant with 
the absolute blue magnitude of the galaxies, we would ex- 
pect a positive correlation between [O n] flux and rest-frame 
blue magnitude. Therefore, if the colours of the [O n] emit- 



ters are constant, we would expect a positive correlation 
between [O n] flux and observed 7^-band magnitude. This 
would make the completeness factor of 0.25 in point (i) an 
underestimate. On the other hand, if the colour term is bluer 
for galaxies with a higher star-formation rate - and hence 
brighter [O n] emission - then this would counter the posi- 
tive correlation. Studies of [O n] emission at z ~ — 1 (e.g. 
Hogg et al. 1998) do not allow us to separate these effects, 
so we adopt the 0.25 completeness factor and discuss the 
effects of extreme departures from this value in Section 3.3. 

3.3 Evidence of clustering 

The distribution of photometric redshifts in the field of 
CVB13 has a clear peak at z = 1.4, as is shown in 
Fig-Ill Spectroscopic follow-up however shows that the single 



10 C. van Breukelen et al. 




Figure 4. The three-dimensional distribution of galaxies within our spectroscopic sample (see Fig. [2] red circles) in comoving coordinates 
(X corresponds to RA, Y to Dec and z to redshift). The plotted box reflects the field of CVB13 as defined in section 2.3; the front and 
back faces are defined by the RA and Dec limits as described in the caption of Fig. [T] The galaxies are colour-coded to redshift, from 
blue at z = 1.25 to red at z = 1.53. Right: XY-projection of the galaxy distribution. 




1.25 1.30 1.35 1.40 1.45 1.50 1.30 1.35 1.40 1.45 1.50 

^spec ^spec 

Figure 5. Left: histogram of the spectroscopic redshifts of the targeted galaxies with detected [O II] emission in the field of CVB13. The 
shaded bin represents the galaxy observed with GMOS. Overplotted is the expected number of [On] emitters per redshift bin, taking 
into account our detection completeness (see section 3.2). The effects of the sky lines are somewhat smoothed in the plot owing to the 
binning of the histogram. Right: The probability of the detected number of [O II] emitters per redshift bin being due to the background 
model, using the method of Brand et al. (2003) described in Section 3.3. The solid black line is for a bin size of Az = 0.0025, the dotted 
blue line for Az = 0.005, and the dashed red line for Az = 0.01. 



redshift-peak at z = 1.4 separates into three less significant 
peaks at z = 1.40, z = 1.45, and z = 1.48, as pictured in 
Fig. [5] (left). There is also a clear overdensity at z = 1.28 
which was not detected in the photometric redshift cata- 
logue. 

Following the method of Brand et al. (2003) we calcu- 
late the probability in each redshift bin (using Poissonian 
low-number statistics) that the number of [O n] emitters 
is greater than the number detected, given the number of 
[O n] emitters predicted by the model. The result is shown 
in Fig. E] (right) for bin sizes of Az = 0.0025 (120 bins, 
black solid line), Az = 0.005 (60 bins, blue dotted line), 
and Az — 0.01 (30 bins, red dashed line). For a bin size 
of Az = 0.01 the redshift peaks at z = 1.28 and z = 1.45 
both have a probability of 10 -5 . As there are 30 bins with 
this bin size, one would expect to find such a low value in 
0.03 per cent drawn from a set of random realizations. Based 
on these statistics we deem the structures at z = 1.28 and 



z — 1.45 to be highly significant overdensities. Other peaks 
in the redshift distribution could reflect genuine structures 
but could alternatively be due to Poisson noise in the back- 
ground distribution. Note that if the completeness factor of 
point (i) in Section 3.2 is smaller than 0.25, the overden- 
sities would be even more significant; if the factor is larger 
than 0.25 this would reduce the probability of the two struc- 
tures up to ~ 0.01 for a factor of 1. However, this large a 
completeness value is disproved by the existence of object 
CVB13_24B which has strong [On] emission but K > 20.6, 
and by the occurence of similar high equivalent width ob- 
jects in the sample of Hogg et al. (1998). 

We do not detect any obvious signs of clustering in RA 
and Dec in our spectroscopic sample (see Fig. |4]). This is 
principally because (i) we were unable to spectroscopically 
observe every photometric candidate cluster galaxy within 
the allocated time (ii) the MOS instrument was unable to 
completely sample the cluster core, as slitlets cannot be 
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Figure 6. Left: The line-of-sight velocity distribution for the structure at z = 1.28. The solid lines in both the histogram and the fitted 
Gaussian comprise all galaxies within a range of ± 1000 km s _1 from the centre. The dashed lines show the sample within ± 2000 km 
s . Bight: Velocity distribution for the structure at z = 1.45. Line styles are as in the diagram on the left. 



placed within 1.5" along the spatial axis, nor next to each 
other along the dispersion axis as the spectra will overlap, 
and (iii) we could only determine redshifts for [O n] emit- 
ting galaxies and not for the red, passively-evolving galax- 
ies, which are customarily more clustered (e.g. Meneux et 
al. 2006). 



3.4 Velocity dispersions and mass estimates 

The line-of-sight velocity distribution of the two statistically 
robust structures at z = 1.28 and z = 1.45 is shown in Fig. [6] 
The bi-weighted mean of the redshifts and the 'gapper' esti- 
mate of the velocity dispersion are calculated for each of the 
structures, using the procedure outlined in Beers, Flynn & 
Gebhardt (1990). Subsequently, we estimate the masses of 
the groups using the relation between cluster mass M200 (the 
mass contained within a sphere of radius r2oo for which the 
mean density is 200 p cr ) and virial velocity found by Evrard 
(2004) (Eq.[5|, substituting ov = V3a v ,i os : 



M200 = 



lO 15 h -1 M ( V3a v 



H/Ho 



( v3ov,iob V 
V 1080 km s" 1 / 



where for a flat Universe: 
H(z) = H y/nA + n M (l + z) 



(2) 



(3) 



An alternative mass estimate is based on the stellar 
light of the groups observed in the if-band. This is done in 
a number of steps. First, we calculate the Tf-band luminosity 
of each galaxy in terms of L* K at the appropriate redshift. 
This parameter is determined by assuming MjJ- = —24.34 for 
the cluster luminosity function (Lin, Mohr & Stanford 2004) 
at z = 0, and passive evolution of the A'-band luminosity 
function given a formation redshift of zf or m = 4. Next, we 
sum all the luminosities of the observed galaxies and cor- 
rect for the missed fraction of light due to the flux limit, 
according to the luminosity function of Lin et al. (2004). To 
arrive at the total mass we then take a mass-to-light ratio 
of (A//M Q )/(L/L ) of 75/i (Rines et al. 2001) which is as- 
sumed constant with redshift in terms of L* K (note that as- 
suming passive evolution of the luminosity function implies 
L* K itself increases with redshift). The mass thus estimated 
is only a lower limit, as we do not correct for the fact that 



our spectroscopic sampling of the group members is most 
likely incomplete. 

As is illustrated in Fig. (left) the mean redshift and 
velocity dispersion of the structure at z = 1.28 is highly de- 
pendent on the limit imposed on the individual galaxy veloc- 
ity for group membership. If we assume all galaxies within a 
range of 2000 kms -1 are members of the group, we obtain 
a mean redshift of z — 1.282 ± 0.002 and a velocity disper- 
sion of <t Vj i os = 668 ± 149 km s" 1 . Placing the limit at 1000 
kms -1 however gives z = 1.2847 ± 0.0008, <r v ,i os = 236 ± 63 
kms - . Assuming the group is virialised, the dynamically 
inferred masses are 5 x 1O 14 M and 3 x 1O 13 M respec- 
tively. The mass estimates from the total stellar light are 
8 x 1O 13 M and 4 x 1O 13 M respectively for the two galaxy 
samples. 

The structure at z — 1.45 gives similar velocity dis- 
tributions for the samples in the two velocity ranges (see 
Fig. E] [right]). We arrive at z = 1.4542 ± 0.0004 and 
o"v,ios = 113 ± 57 kms -1 ; the inferred dynamical mass is 
only 5 x 1O 12 M . The mass derived from the stellar light 
yields a much higher estimate of 9 x 1O 13 M . Since we have 
spectroscopic confirmation for six to eight galaxies in the 
group with a total stellar mass of 1.6 — 1.8 x 1O 12 M we 
believe ~ a few 1O 13 M is a sensible minimum mass for the 
associated dark matter halo. 



3.5 X-ray and radio properties 

We have performed an analysis of the X-ray emission around 
CVB13. Fig. [7] (left) shows a three-colour (Bi'K) image 
of the central 3' x 3' field with X-ray and radio (VLA B- 
array) contours overlaid. We detect X-ray emission at RA = 
02 ,l 18 m 09 a and Dec = -05°00'30" (object X2, see Tabled 
and Ueda et al. [in prep] source no. 0712) which is coin- 
cident with the optical/near-infrared cluster position. The 
X-ray source appears extended with a deconvolved size of 
~ 5" (the XMM point spread function is ~ 6"). The total 
X-ray flux is Fx = 9.1 x 10 -18 Wm -2 , which would mean 
a luminosity of Lx = 8.9 x 10 36 W if the emission is asso- 
ciated with the cluster at z = 1.28 or L x = 1.2 x 10 37 W 
at z = 1.45. This would be indicative of a cluster mass of 
~ 1O 14 M for one of the clusters, or ~ 8 x 1O 13 M each if 
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Figure 7. Left: Bi' K colour image (Furusawa et al. [in prep.]) of the central 3' X 3' region of CVB13. X-ray contours are from a 
signal-to-noise radio map in the 0.5-2 keV band and are overlaid in purple (contour levels at [\/2(T, 2a, 2\/2<j, ...]); the 1.4-GHz VLA 
B-array radio map is contoured in blue (contour levels at [4<r, 4\/2cr, 8<r, ...]). The radio beam size is 5" X 4" at PA = 170°. X-ray and 
radio sources are labelled as in TablefS] Right: Colour image as left, zoomed in on the central 1' X 1' region. The radio contours are from 
the higher resolution VLA A-array radio map (contour levels at [2%/2ct, 4<t, 4\/2o", •••]), with a beam size of 1.9" by 1.6" at PA = 22°. 



ID RA Dec. F z vhoto z spcc Ref. 

[h m s] [° ' "] [W m" 2 ] 



XI 02:18:06.52 -05:00:43.9 2.5 ± 0.3 X 10~ 17 1.06 

X2 02:18:08.37 -05:00:30.3 9.1 ± 2.1 X 10~ 18 

X3 02:18:13.83 -04:59:39.2 6.4 ± 1.5 X 10" 18 - 1.258 Simpson et al. (in prep) 



[h m s] [° ' "] [/.Jy] 



RBI 


02:18:03.66 


-04:59:30.1 


122 ± 26 


1.05 






RB2* 


02:18:06.83 


-04:59:16.8 


174 ± 50 


1.11 


1.276 


This paper (CVB13.5) 


RB3 


02:18:07.94 


-05:01:45.0 


256 ± 27 








RB4 


02:18:08.59 


-05:00:40.0 


102 ± 29 




0.493 


Simpson (priv. comm.) 


RB5 


02:18:09.48 


-04:59:45.5 


2347 ± 28 


1.35 


1.094 


Simpson et al. (in prep) 


RB6 


02:18:09.55 


-05:02:00.0 


68 ± 22 




1.268 


Simpson (priv. comm.) 


RB7 


02:18:11.89 


-05:01:29.6 


49 ± 21 


0.75 


0.918 


This paper 


RB8 


02:18:12.44 


-05:02:01.3 


104 ± 27 


1.32 


1.453 


This paper (CVB13.19) 


RA1 


02:18:08.30 


-05:00:27.2 


67 ± 20 


1.02 






RA2 


02:18:08.57 


-05:00:52.8 


73 ± 30 


0.94 






RA3 


02:18:09.10 


-05:00:27.2 


66 ± 29 


1.27 


1.454 


This paper (CVB13.21) 


RA4 


02:18:09.23 


-05:00:30.0 


109 ± 38 


1.32 







Table 2. Table of the X-ray and radio sources in the 3' X 3' central region of CVB13. The IDs correspond to the labels in Fig. [7] The RA 
and Dec are given in columns 2 and 3. The values in column 4 are the X-ray flux (0.5 - 10 keV) for the X-ray sources, and the 1.4-GHz 
radio flux density for the radio sources. The photometric redshifts (column 5) were determined as described in VB06. The spectroscopic 
redshifts (column 6) are from various sources as referenced in column 7. *This object appears to be a resolved source in the B-array data, 
but the A-array data shows that it actually consists of 3 separate objects. The central object on its own has a flux of F = 113 ± 24/xJy. 



both clusters contribute equally to the X-ray flux (Reiprich 
& Bohringer 2002). The signal-to-noise ratio is too low to 
allow us to detect the 6.4 keV iron line and thus we can- 
not assign a secure redshift to the X-ray source. However, 
enough photons are detected in the source to conduct a 
crude hardness ratio (HR) spectral analysis. The HR-values 



are HR2 = (Gs - C 2 )/{C 2 + G) = 0.38 ± 0.22 and HR3 
= (C 4 - G)/(G + G) = 0.01 ± 0.16 where G is the 0.5 
- 2.0 keV count rate, G the 2.0 - 4.5 keV count rate, and 
C4 the 4.5 - 10.0 keV count rate (Ueda et al. [in prep.]). 
These ratios mean the X-ray spectrum is much harder than 
typical for a cluster, and is more characteristic for an ab- 
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sorbed active galactic nucleus (e.g. Delia Ceca et al. 2004). 
However, we do not find evidence for a single obscured AGN 
(e.g. a 24-/Ltm or radio source) at the exact position of the 
X-ray emission. This is a puzzling system which is at face 
value an extended hard X-ray source. It requires an emission 
mechanism beyond any simple addition of a single obscured, 
non-thermal AGN and thermal cluster emission. Two such 
mechanisms that can cause the extended hard X-ray emis- 
sion are inverse Compton scattering (see for example Fabian 
et al. 2006b) or confusion of more than one hard-spectrum 
X-ray source. 

There is a significantly brighter X-ray source ~ 40 arc- 
sec to the southwest of CVB13 (XI, see Fig. \7\ [left]). It 
is uncertain if it is physically associated with the weaker 
extended X-ray emission; however it is certainly associated 
with a distant galaxy. A third X-ray source in the field (X3) 
is associated with an AGN at a redshift coincident with the 
2 = 1.28 structure (see Table [2}. Table [2] also shows the 
properties of the eight brightest radio sources in the field 
shown in Fig. [7] (left). There are radio sources associated 
with both the z = 1.28 and the z = 1.45 clusters, but the 
bright radio source RB4 closest to CVB13 is associated with 
a foreground galaxy at z ~ 0.5. In Fig. [7] (right) we plot a 
zoomed-in view of the CVB13 field and overplot the con- 
tours of the VLA A-configuration map. This confirms four 
weak radio sources tentatively seen in the B-configuration 
map (Fig. [7J [left]). One of these objects (RA3) has spectro- 
scopically been confirmed to be at z = 1.45. RA4 has similar 
colours and photometric redshift and is therefore likely to 
be part of the same cluster. There is inconclusive evidence 
that the radio source associated with this galaxy has the 
'head-tail' structure characteristic of galaxies moving with 
respect to an intracluster medium. 

If the radio emission is solely due to synchrotron ra- 
diation from a starburst, the SFR of a radio source with 
a 1.4-GHz flux density of 100 ^iJy is ~ 14OOM yr _1 in 
stars with masses of > 5Mq (see Condon 1992). Assum- 
ing a Salpeter IMF, this would mean a SFR(M > 0.1 M ) 
of ~ 7000 M yr" 1 . The SFR derived from the [O n] emis- 
sion of these objects is only ~ 10 M yr -1 , and although this 
estimate does not take extinction into account and is there- 
fore a lower limit, the sub-millimetre SCUBA HAlf Degree 
Extragalactic Survey (SHADES, e.g. Coppin et al. 2006) 
covers this region and does not identify any highly obscured 
starbursts. If, on the other hand, all radio emission is as- 
cribed to AGN activity we would expect (using Simpson et 
al. 2006, Eq. 9) a total X-ray luminosity of ~ 10" 16 Wm" 2 
in the 2-10 keV range. The object X2 only has a flux in this 
band of ~ 10~ 17 Wm -2 ; therefore confusion between hard 
X-ray sources associated with the radio sources is a possible 
explanation for this seemingly extended hard X-ray source. 



4 CONCLUDING REMARKS 

We have obtained multi-object spectroscopy on a photomet- 
rically selected cluster candidate at z — 1.4. Instead of a 
single massive cluster, we find three projected structures of 
galaxies at z = 1.40, z = 1.45, and z = 1.48, one of which 
(at z — 1.45) is statistically robust. We also serendipitously 
detect a robust structure at z = 1.28. Both the structure at 



z = 1.45 and the one at z = 1.28 have masses ~ 10 M 
and may therefore be termed poor galaxy clusters. 

The evolutionary status of both the z — 1.28 and 
z = 1.45 clusters is uncertain. The cluster at z = 1.45 has 
a significantly lower virial mass than its mass derived from 
the total observed luminosity. It is very likely that this struc- 
ture is not yet virialised, which would suppress the derived 
(three-dimensional) velocity dispersion of ~ 200kms _1 . 
The velocity shear (half-width to zero intensity) observed 
in CVB13_21 - located centrally within the cluster - is 
320 km s -1 , which shows that the measured velocity disper- 
sion of the cluster is unlikely to represent the total cluster 
mass. Note however that low-number statistics make an ac- 
curate determination of the velocity dispersion unfeasible. 
The total measured SFRs are ~ 45 and ~ 6OM yr _1 re- 
spectively for the clusters at z — 1.28 and z = 1.45, which 
for a mass of ~ 1 x 1O 14 M within 4x4x4 Mpc 3 (comoving) 
is a factor of 3 - 5 higher than typical star-formation rates 
observed in the field at these redshifts (e.g. Conolly et al. 
1997) 

The analysis of the cluster galaxies in this paper is based 
solely on redshifts estimated from [On] emission line fea- 
tures. Therefore the overdensities we find at z = 1.28 and 
z = 1.45 only consist of star- forming galaxies, whereas the 
candidate cluster found by VB06 was selected with a photo- 
metric redshift algorithm which is sensitive to red, passively- 
evolving galaxies. Of the 18 highest-priority targets of sam- 
ple 1, 10 were observed, and only six redshifts were obtained, 
two of which were confirmed to be cluster members. The 
other galaxies were all priority 3 targets of sample 1. We 
suspect there will be an appreciable number of red galax- 
ies in the cluster that we were unable to confirm spectro- 
scopically. More observations of the central cluster galaxies, 
specifically for the cluster at z = 1.45, may also produce a 
higher velocity dispersion measurement. To thoroughly un- 
derstand the nature and masses of the structures, we need 
to be able to obtain redshifts for this type of galaxy as well. 
Deep, near-infrared spectroscopy, for example with the Fi- 
bre Multi Object Spectrograph (FMOS, Dalton et al. 2006), 
will give a more complete sample and therefore yield a more 
comprehensive result for high-redshift galaxy clusters. 

The extended X-ray luminosity at the position of these 
structures is characteristic of a cluster with a mass of ~ 
1O 14 M . However the X-ray spectrum is too hard to simply 
be caused by thermal cluster gas emission. A non-thermal 
effect like inverse Compton scattering could cause an X-ray 
detection with a spectral shape reflecting the energy dis- 
tribution of the underlying electron population. The most 
likely explanation, however, is confusion caused by more 
than one AGN contributing to the X-ray emission. To shed 
light on the exact nature of this system, it will be interest- 
ing to study the cluster gas via the Sunyaev-Zel'dovich effect 
and the dark-matter mass via lensing. 

The original mass estimate of CVB13 by VB06 of 
3 x 1O 14 M is likely to be overestimated due to the su- 
perposition of the two lower-mass structures at z = 1.40 
and z — 1.48. The single structure at z = 1.28 contains 
too few galaxies to be isolated by the same selection algo- 
rithm. This result suggests that spectroscopic follow-up is 
a vital element of photometric cluster surveys as the latter 
are prone to projection effects and can seriously overestimate 
cluster masses. To determine the impact of these effects, we 
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calculate the probability of finding three groups in projec- 
tion within Az — 0.1 in a survey such as the one presented 
in VB06. Assuming a Sheth-Tormen cluster mass function 
(Sheth & Tormen 1999) we expect to find 121 groups with 
masses of a few 1O 13 M0 and 6 clusters with masses of 
~ 1O 14 M0 within the redshift range 1.25 ^ z ^ 1.55 in a 
field of 0.5 deg 2 . Neglecting spatial clustering of the groups, 
there are on average two superpositions of three groups with 
masses ~ 1O 13 M within a box of 2 x 2 Mpc 2 and Az = 0.1. 
The probability of finding a projection of two groups and one 
cluster of ~ 1O 14 M (as for CVB13) is 45 per cent. Running 
the cluster algorithm of VB06 on a simulated superposition 
of three groups shows that in order to distinguish between 
massive clusters at this redshift and projected groups with 
a redshift separation of Az — 0.05, the photometric redshift 
error needs to be as small as = 0.005. The COMBO- 

l + z 

17 survey (Wolf et al. 2004), which determined photometric 
redshifts from photometry in 17 bands, obtained a best red- 
shift error of = 0.01 for galaxies with R < 21 (c.f. 
R ~ 24 for CVB13 targets). Thus attaining the order of 
accuracy needed in photometric redshifts to determine the 
cluster mass function at high redshifts seems very challeng- 
ing, and we caution against attempting this with photomet- 
ric techniques alone. 
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